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We investigate the phenomenon of fermionic pairing with mismatched Fermi sur- 
faces in a two-species system in the presence of Feshbach resonance, where the 
resonantly-paired fermions combine to form bosonic molecules. We observe that 
the Feshbach parameters control the critical temperature of the gapped BCS super- 
fluid state, and also determine the range over which a gapless breached pair state 
may exist. Demanding the positivity of the superfluid density, it is shown that al- 
though a breached pair state with two Fermi surfaces is always unstable, its single 
Fermi-surface counterpart can be stable if the chemical potentials of the two pairing 
species have opposite signs. This condition is satisfied only over a narrow region in 
the BEC side, characterized by an upper and a lower limit for the magnetic field. 
We estimate these limits for a mixture of two hyperfine states of 6 Li using recent 
experimental data. 

PACS numbers: 03.75. Ss, 74.20.-z, 05.30.Fk, 03.75.Kk 



1. INTRODUCTION 

The problem of fermionic pairing and superfluidity with mismatched Fermi surfaces has 
been widely investigated in recent years from theoretical standpoints 



9. 



10 



11 



12 



13] as well as experimentally [33j, |34|, [35] 



A mismatch in Fermi surfaces can be easily realized in a two-species fermionic system 
where the pairing species have unequal populations or different masses/ chemical potentials. 
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Typically, such a system would consist of two different fermionic atoms (e.g., 6 Li and 40 K) 
or alternatively, two hyperfme states of the same atom (e.g., states \F = 1/2, mp = 1/2) 
and \F = l/2,m F = -1/2) of 6 Li atoms). 

Several phases have been proposed to describe the possible ground state of such a system, 



including the FFLO phase jl6, [l?], 



2(]] (pairing with non-zero centre-of-mass momentum, 



where the order parameter shows a spatial variation), the gapless BP (breached pair) phase 
(pairing with zero centre-of mass momentum: where the order parameter is non 
zero but the excitation energy becomes zero), and the inhomogeneous phase-separated state 
0,S], where any two pure states coexist. 

In this note, our focus is on the homogeneous and gapless breached pair state, also 
known as the Sarma state. Sarma [15], in the early studies of superconductivity, predicted 
a spatially isotropic, homogenous and uniform state with gapless excitation modes in the 
presence of a magnetic field. However, for weak coupling BCS theory, this gapless breached 
pair state marks the maximum of the thermodynamic potential, and thus, cannot be the 
stable ground state of the system. This is the well-known Sarma instability. In the last few 
years, several mechanisms were put forward to avoid this instability. According to Forbes et 



al. 18[, a stable Sarma state is possible in a model with finite range interaction where the 
momentum dependence of the pairing gap cures the instability. It has also been proposed 
by a number of workers jf], Q, 20, 21] that the breached pair state becomes stable in the 



deep BEC regime, if the BCS-BEC crossover picture is taken into account. He et al.,in a 
very recent work [23|] has argued that the breached pair state can be a possible ground state 
in the weak coupling region for a two-band Fermi system. 

To study the breached pair state, we start with a two-species fermionic system. In 
addition to the weak BCS attraction ( denoted by —gi), we consider a strong interaction 
((72) of the Feshbach variety which couples a fermion of type a with a b fermion to form a 



bosonic molecule B. Our model resembles the one used in 



24 



25j , but we extend it to cover 



the two-species case. The system is described by the Hamiltonian: 

h = J2( 2 "-^b)bIb +j2 ^4°p+Z) J2 a l> b - P > b -P a P+92 J2i B l a P b - P + a l b -P B o] 

(i) 

Here a p , a p are the creation and annihilation operators for atom a, while b p , bp are the 
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corresponding operators for atom b. Also, e p = e p — fi a and e b = e b — while e p , e b p are the 
respective kinetic energies and /i a , ^ the respective chemical potentials for species a and 
b. The annihilation and creation operators for the composite boson B are B and (We 
restrict ourselves to the case where only zero- momentum bosons are formed). The chemical 
potential for the bosons is fiB, and 2z/ is the threshold energy of the composite bose particle 
energy band. 



Using this Hamiltonian, we study the effect of the Feshbach parameters first on the gapped 
BCS superfluid state, and then on the gapless Sarma state, taking a variational approach. 
The stabilit y o f such a gapless state is analysed ensuring the positivity of the superfluid 
density 0, 29]. Considering the BCS-BEC crossover picture, we show that the breached 
pair state is stable only in a narrow region in the BEC side, bounded by two magnetic field 
values. 



2. GAP EQUATION AND GAPLESS EXCITATIONS 

We want to study the ground state of the two-fermion system by a variational method. 
We take |\&) = \F) ® \B) as the trial form of the ground state of the system. Here \F) 
= \BCS) =n{U p + V p dU)_ p )\0) , i.e, the probability of the pair (a p f, 6_ p |)being occupied is 
|V^,| 2 , and the probability that it is unoccupied is \U P \ 2 = 1 — |V^| 2 , while \B) is the ground 
state for the condensate part of the boson subsystem. Our variational prescription is that 
\B) has to be chosen in such a manner that it is an eigenstate of the annihilation operator 
B, which would make it a coherent state. In terms of Fock states \B) =HC n \n). Using the 
normalization condition, we arrive at \B) = exp(— a 2 /2)exp(a5^)|0). The ground state of 
the system is consequently given by 

l*> = Y[( u p + V P a l b -p)\°) ® exp(-« 2 /2 + aB^)\0) (2) 

where a = V^b, ^Yb being the expectation value of the total number of bosons in the 
condensed state. 

From the Hamiltonian (TfJ and the ground state wave function (j2J), the ground state 
energy of the system would be 
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E = 5>1 + ~4K 2 ~9iJ2 U p V p U p' V p' + ( 2 » - M« 2 + 2^ 2 a Yl U p V p ( 3 ) 



Minimizing with respect to V p and a we get 

4ejy p - 2g eff £ 0p^(l7 p - = (4) 

p,p' 

where e+ = (e£ + eJ)/2. 

(7 e // is defined by the relation 

g e ff = gi + gi/(2v-LiB) (5) 



We note that this expression matches with the one obtained by Ohashi et al [24( using 
diagrammatics. In order to avoid a possible ultraviolet divergence, the summation in the 
right hand side of equation(j4]) has to be carried over upto a cutoff. For a metallic super- 
conductor, this cutoff is played by hun, being the Debye frequency In ultracold atomic 
systems, the cutoff is determined by the range of the interatomic potential. 

Now, if we choose g e ff J2 p > U p iV p i = A, the usual form of the gap equation follows, and 
we can identify A as the gap in the excitation spectrum. A similar relation is obtained 
in finite temperature systems as well, if we incorporate the appropriate Fermi distribution 
functions in the expression of the free energy. In terms of this gap parameter, U p and V p 
can be expressed as ift = \ [ 1 H — , = | and V% = | ( 1 , e l ) . 

The critical temperature Tq (temperature at which the gap vanishes) is proportional to 
A , the BCS gap in the weak coupling limit, which is given by 27] 
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A = — e F exp( — ) (6) 

e 2 P{0)9eff 

Here p(0) is the density of states in the Fermi level, e is the base of natural logarithms, and 
ep is the Fermi energy. Thus the value of A ( and T c ) can be raised or lowered by adjusting 
2u, the tunable parameter depending on the Feshbach resonance process. 



Next we study gapless excitations 
BCS-like treatment is of the form 



he quasiparticle dispersions as obtained from standard 



14j: 



„ 4 n - et I r» - c>; 

"p 



Ef = ^ + V (^T^) 2 + A2 ( 7 ) 
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FIG. 1: (E%> -p) curve when the Feshbach term is absent 
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FIG. 2: (Ep ,b -p) curve for — 2u = .5, §2 = .2 
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FIG. 3: (Ep b -p) curve for \Ib — 2v = —.5, §2 = .2 



Depending on the values of particle masses and corresponding chemical potentials, these 
quasiparticle excitations can be negative, thus leading to gapless excitations. This is possible 
only if the magnitude of A is less than a critical value A c 3l. Il4|. 



\m b fx b - m a jj a \ 
2y/m a m b 



(8) 



When I A I > A c , both E% and E h remain positive for all values of p. This corresponds to 



usual BCS pairing. When |A| < A c , either E't or E b crosses zero at the points: 
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Pi,2 = ( m bHb + mafia) T [{m h fi h - m a /i a ) 2 - 4m a m b A 2 ] 1/2 (9) 

The difference between pi and j>2 gives the span over which we get a gapless region in the 
parameter space. The state with gapless excitations marks the coexistence of the superfluid 
and normal components at zero temperature, and is called the Sarma phase, or the Breached 
Pair state. When we include the Feshbach term in the Hamiltonian, we have a control over 
this Sarma phase as well. 

In figures(l,2,3), we plot the E-p curve for a two-species system. Here we scale all 
energies by e bp ( Fermi energy of species b) and all momenta by p bp (Fermi momentum of 
species h). It is evident that in this convention, m b = .5 and fib = l(in the BCS limit). We 
choose m a = .1, fi a — 6 and multiply quantities gi, g 2 , (2z/ — fi B ) and g e ff by p(0), the 
density of states at the Fermi level to get dimensionless quantities g~\, g~2, (2z/ — fis) and 
geff- Let g~i = 0.3. Had there been no Feshbach coupling g 2 , we would get a gapless region 
from p=1.01 to p=1.85 ( in units of pb F ) as seen from Figure-1. 

If 9~2 — -2 is introduced in the system, and we choose v in such a way that — 2z/) = 0.5, 
we get the Sarma phase for a wider region, from p=1.00 to p=1.95 (in units of pt F ) as seen 
from Figure- 2. 

If, on the other hand, (/i~b — 2z/) = —.5, the gapless phase vanishes entirely (Figure-3). 



3. STABILITY ANALYSIS 



The stability of a superfluid phase has been studied in different ways like demanding 



21 



ZM , ensuring non-negative eigenvalues of the number 



3l| , minimization of thermodynamic potential 



positivity of the superfluid density 
susceptibility matrix of the system 
etc. 

Here we adopt the first criterion and demand that the superfluid density n s must be 
positive. If the fermions are charged, this is equivalent to saying that the Meissner mass- 
squared must be positive [3l|, since n s , the superfluid density, and M 2 , the Meissner mass 



squared obey the relation : M 2 



q being the electronic charge and m, the mass 



of the particle. In our treatment, we work with the superfluid density directly, since we 
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are talking of a charge neutral Fermi system. However, the term "Meissner mass" can still 
be used just to continue the analogy with superconductors, as it only serves to express the 
superfluid density to within a multiplicative constant. 

For convenience, we assume that the two species have equal masses. Therefore, equation 
(Q takes a simpler form: 

p 1>2 = 2m[fL ± W- A 2 ] 1 (10) 
where /2 = + £fo)/2 and 5/i = (/i a — /Ub)/2. Thus, A c = 5fi here. 



He et al 



28 



29| have shown that when the two species have equal masses, the superfluid 



density can be expressed as 



;/ __ m r^^^-A) (n) 



V^/i 2 - A 2 



p'I + pI 



where r\ = - — . We note that here n marks the bare fermion density or the density of 

<oix 1 n 

the atoms that have not been part of the condensate yet ( i.e, unpaired fermions, fermions 

forming the Cooper pair and fermions that constitute the non-condensate bosons all are 

counted in n). Since the total number of bare Fermi atoms is conserved, so n + 2Nb= 

constant, Nb being the expectation value of the total number of bosons in the condensate. 

If A < A c , for n s is to be positive, — = has to be less than 1, which is satisfied if 

V V - A 2 

A 2 

rj < 1, and 77 — 1 > -r— 

d/r 

Now when both pi and p 2 are real and there is a breached pair phase between them (this 
state has sometimes been termed as BP2 state in the literature 37J, [38]), the bare fermion 
density can be written as: 

1 P 2 o , 1 



n = 7T~2 / P 2 dp+^ 
2tt z L, tt z 



rpi roo 

/ v p 2 P 2 dp+ / v p Yd P 

JO Jvo 



(12) 

pi " L^o Jp2 

Here the first term denotes the contribution from the gapless region between p\ and p 2 , 

i.e, the normal component. The next two integrals take care of the contributions from the 

gapped superfluid regimes, one from momenta to p\ and the other from p 2 to infinity. The 

terms, when rearranged, gives 

pI+pI i 1 , r\ 2TT2 r °° 

67T 2 7T 2 

In the weak coupling limit, the last two integrals are very small, and hence, can be 
neglected. Therefore, 77=1 , and n s is negative. Thus, the Sarma state is unstable here. 



+ ^(-/ P U P d P + J PVpdp) (13) 



8 




We shift our focus to a special case, where only one of p\ and p<i is real ( in the lit- 



erature, this is called the BP1 state 
and the bare fermion density is 



33, 



38|). In this case the E-p curve resembles figure 4, 



n 



I '-2 ' /' -2T.-2 



6lT 2 7T Z 



+ - 2 p\ 2 d P (14) 



This yields r\ < 1, and thus the stability criterion is fulfilled. Now, from equation ([10]), p\ is 
imaginary if a/ <5/i 2 — A 2 > /2, i.e, — /i a /i{, > A 2 . So, Sarma phase is stable only in a region 
where the chemical potential of one species is positive, and the other, negative, provided 
the magnitude of their product is greater than A 2 . 



Leggett has shown [32j that in the BCS-BEC crossover picture, the chemical potential 
can be determined by solving the gap and the number equations. Extending these equations 
to the two species case, we get 

°° 1 1 m 

^7 P ~ y /^+W ) = 2^a s (15a) 
^ Jt 2 + A 2j 3vr 2 1 ; 

P=P2 V P 

°° en P2 en k 3 

where a s is the a — b scattering length and k ap and kb F correspond to the Fermi wave 
number of the more and the less populated species respectively. Let k' denote the wave 
number corresponding to the breaching point for the BP1 state, i.e, P2 = hk'. Converting 
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H > 



BCS side 



FIG. 5: Behavior of (i a and fib in the crossover picture 



the sums into integrals, we obtain that in the weak coupling limit, fi = (e ap + e& F )/2 
h 2 (k 2 F + kf )/Am, as expected. In contrast, in the strong coupling limit, 

h 2 2e bF (k bF a s 



/i 



1 + 



2k' 3 a 



(16) 



2ma 2 3tt \ n 

provided k' 2 « p, (a condition which is satisfied if the population imbalance is small 
compared to the total population). In this case p, asymptotically approaches —h 2 /2ma 2 , i.e, 
half the binding energy of the molecule. 

If we think of \i a and /zj, separately, the first one differs from the other by the Fermi energy of 
the excess fermions. Moreover, in the presence of a magnetic field H, there is an asymmetry 
between the chemical potentials given by m^if, m# being the fermion magneton. So we 
have 

f " = K^ +m ' J/ ) (17) 
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Therefore, 



ft 2 2e bF (k bp a s 



+ 



1 



7T 



± 



h 2 k' 2 m B H\ 

-a — + — 

Am 



) 



18) 



2ma 2 ' 37r 

Our domain of interest is when ii a is positive, and fi b is negative, i.e, while approaching 
the BEC side, \i b has already crossed zero but /i a has not. Now, \i a and \i b becomes zero at 
magnetic field values H% and H 2 respectively, as sketched in figure (5), where 



H x = (2/m B ) 
H 2 = (2/m B 



h 2 



si , 



2ma 2 sl 



h 2 



2ma 2 s2 



+ 



3tt 

^b F {kb F 0's2) 



3tt 



1 





h 2 k 12 ' 




Am 


2k^ s2 \ 


h 2 k' 2 ' 




Am 



(19a) 
(19b) 



Here a s i, a s2 are the respective values of the scattering length at Hi and H 2 . Between these 
two magnetic field values, the Sarma state will be stable. 

4. ESTIMATES OF Hi AND H 2 



Experiments with population-imbalanced fermionic systems have been done by Zwierlein 
et al. [33], 34] and Partridge et al 35]. They obtained the signature of superfluidity in an 
unequal mixture of two spin states of 6 Li atoms, and a quantum phase transition between 
the superfluid state and the normal state was observed at a critical polarization. Another 
method for experimental detection of the breached pair phase has been suggested by Yi et 



al. 



371 ] . However, no clear signature of this gapless phase has been obtained till date. 
We now use the data obtained from these experiments on population-imbalanced gas of 6 Li 



atoms 



33 



34j to make an estimate of the magnetic field values corresponding to the breached 



pair Sarma state. To be able to use the expressions fll9alll9bl) . we need to know the scattering 
length as a function of the magnetic field. This is provided by a s = clq \ 1 — — — 1 ■ whoro 



H-H 0/ 

ao is the background scattering length, T is the width of the resonance, and H is the position 
of the resonance peak. 

We use this expression for scattering length, and put a Q = 45.5r (r =Bohr radius), 
which is the singlet scattering length for 6 Li. We also take n a = 1.8 x 10 7 , rib = 2.6 x 
10 6 as population of the two species, H = 834 G, T= 300G ( These values correspond to 



experimental data for Feshbach resonances in 6 Li as reported in 13c 



34 



36j.) We now solve 
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equations (I19a019bl) explicitly to find Hi to be 832.40 G and H 2 to be 833.95 G. It would 
be interesting to speculate that the superfluid observed by Zwierlein et al. in this range is 
of the BP1 variety. 

This estimate of Hi and H 2 is not a highly accurate one. In fact, the range of the breached 
pair state should get shifted a bit towards lower values of the magnetic field. Actually, for 
deriving equations ( I19alll9bl) . we had assumed the magnitudes of the chemical potentials to 
be very large compared to A and kf, an assumption which does not hold at points where 
fi a and fib become zero. However, the values are not unreasonable in view of the fact that 
in the BCS-BEC crossover picture, the chemical potential falls quite sharply after crossing 
zero and quickly becomes a large negative quantity. 



5. SUMMARY AND DISCUSSION 



Here we have studied a two-species fermionic system in the presence of Feshbach reso- 
nance, taking a variational route with an explicit construct of the ground state. The gapless 
breached pair state was discussed, and its stability was analysed ensuring the positivity 
of the superfluid density. We showed that a breached pair state with two Fermi surfaces 
is always unstable, while its single Fermi-surface counterpart is stable when the chemical 
potential of the two pairing species bear opposite signs. 

The stability of the breached pair state is indeed a widely debated issue. Although it 
has often been suggested that the BP1 state might be stable in deep BEC region, nothing, 
to the best of our knowledge, was said anything about how 'deep' that really is. In this 
note, however, we observe that, the requirement that \i a and fib should be of opposite signs, 
automatically puts two bounds in terms of the Feshbach magnetic fields, between which the 
gapless state is stable. Moreover, this stable breached pair state is obtained not in deep 
BEC, but near the vicinity of the point when the average chemical potential crosses zero, 
i.e, right after the onset of condensation. 

Gubankova et al {mJ , while discussing the stability of breached pair states by analysing 
the number susceptibility, commented that stable gapless states with a single Fermi surface 



exist for negative average chemical potential. In a recent paper, A. Mishra et al. 



221 ] reached 



the same conclusion by comparing the thermodynamic potentials of the condensed phase 
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and the normal phase. Although the criterion they arrive at (the negativity of the average 
chemical potential) does not fully match with ours ( the chemical potential of the two species 
to have opposite signs), there is definitely a region of overlap. 

In their experiment, Zwierlein et al. observed that superfluidity breaks down when the 
pairing gap A becomes small compared to the chemical potential difference \i a — fif,. We 
note that this matches with the stability criterion for the BP2 state, since, if A < 8{i and 
T) — 1, the state becomes unstable, as seen from equation ffTTl) . 

As for the BP1 state, we have shown that this state is stable in a region where /i a and 
/ib have opposite signs, which can be achieved by keeping the system between two specific 
magnetic field values. For simplicity we took the two species to have same masses in 
our calculation, but the treatment should be extendable to a situation where the fermion 
species are of different masses, as for example in a 6 Li- 40 K system. 



ACKNOWLEDGEMENT 

The author would like to express her gratitude to her supervisor Prof. J. K. Bhattacharjee 
for his keen interest and helpful comments during the course of this work. Constructive 
suggestions from the referee are also thankfully acknowledged. 



[1] W. V. Liu and F. Wilczek, Phys. Rev. Lett. 90, 047002 (2003). 

[2] P. F. Bedaque, H. Caldas, and G. Rupak, Phys. Rev. Lett. 91, 247002 (2003). 

[3] S-T. Wu and S. Yip, Phys. Rev. A 67, 053603 (2003) 

[4] H. Caldas, C.W. Morais and A.L. Mota, Phys. Rev. D 72, 045008 (2005) 

[5] J. Carlson and S. Reddy, Phys. Rev. Lett. 95, 060401 (2005). 

[6] D. T. Son and M. A. Stephanov, Phys. Rev. A 74 013614 (2006). 

[7] D. E. Sheehy and L. Radzihovsky, Phys. Rev. Lett. 96, 060401 (2006). 

[8] D. E. Sheehy and L. Radzihovsky, Annals of Physics 322, 1790 (2007) 

[9] C. C. Chien, Q. Chen, Y. He and K. Levin, Phys Rev Lett. 98, 110404, (2007) 



13 



[10] Q. Chen, Y. He, C. C. Chien and K. Levin, Phys Rev B 75, 014521, (2007) 
[11] T. Koponen, J. Kinnunen, J. P. Martikainen, L.M. Jensen, P. Torma, New Journal of Physics 
B 8, 179 (2006) 

[12] M. Iskin and C. A. R. Sa de Melo, Phys. Rev. Lett. 97, 100404 (2006) 

[13] K. B. Gubbels, M. W. J. Romans, and H. T. C. Stoof, Phys. Rev. Lett. 97, 210402 (2006) 

[14] Pairing in Fermionic Systems, edited by A. Sedrakian, J. W. Clark, and M. Alford, World 

Scientific Publishing Co, Singapore, 2006 
[15] G.Sarma, J.Phys.Chem.Solid 24,1029(1963) 
[16] A. I. Larkin, Yu. N. Ovchinnikov, Sov. Phys. JETP 20, (1965) 
[17] P. Fulde, R. A. Ferrel, Phys. Rev. A, 135, 550, (1964) 

[18] M. M. Forbes, E. Gubankova, W. V. Liu and F. Wilczek, Phys. Rev. Lett. 94, 017001(2005) 
[19] H.Hu and X.J. Liu, Phys. Rev. A ,73, 051603(R)(2006) 

[20] M. Mannarelli, G. Nardulli and M. Ruggieri, Phys. Rev. A.74, 033606(2006); 

[21] C.-H. Pao, S.-T. Wu, and S.-K. Yip, Phys. Rev. B 73, 132506 (2006). 

[22] A. Mishra and H. Mishra The Eur. Phys. J. D., 53,75,(2009) 

[23] L. He and P. Zhuang , Phys. Rev. B 79, 024511 (2009) 

[24] Y. Ohashi and A. Griffin, Phys. Rev. Lett. 89, 130402 (2002) 

[25] M. Holland, S.J.J.M.F. Kokkelmans, M.L. Chiofalo, R. Walser, Phys. Rev. Lett.87, 120406 
(2001) 

[26] R.P. Feynman, Statistical Mechanics, A Set of Lectures, Addison- Wesley, New York, 1988 
[27] C. J. Pethick and H. Smith, Bose-Einstein Condensation in Dilute Gases, Cambridge Univ. 
Press 2002 

[28] L. He, M. Jin, P. Zhuang, Phys. Rev. B 73, 214527 (2006) 

[29] L. He, M. Jin, P. Zhuang, Phys. Rev. B 73,024511 (2006) 

[30] M. Kitazawa, D.H. Rischke, I. Shovkovy, Phys. Lett. B 637, 367 (2006) 

[31] E. Gubankova, A. Schmitt, F. Wilczek, Phys. Rev. B 74, 064505 (2006) 

[32] A.J. Leggett, Quantum Liquids, Bose Condensation and Cooper pairing in condensed-matter 

systems, Oxford University Press, Oxford, 2006 
[33] M. W. Zwierlein, A. Schirotzek, C. H. Schunck and W. Ketterle, Science 311, 492 (2006) 
[34] M. W. Zwierlein, C.H. Schunck, A. Schirotzek, and W. Ketterle, Nature( London) 442, 54 



14 



(2006) 

[35] G.B. Partridge, W. Li, R. I. Kamar, Y. A. Liao, and R. G. Hulet, Science 311, 503 (2006) 
[36] C. H. Schunck, M. W. Zwierlein, C. A. Stan, S. M. F. Raupach, and W. Ketterle, A. Simoni, 

E. Tiesinga, C. J. Williams, P. S. Julienne, PRA 71, 045601(2005) 
[37] W. Yi and L.-M. Duan, Phys. Rev. Lett. 97, 120401 (2006) 
[38] H. Caldas and A.L. Mota, J Stat. Mech, P08013 (2008) 



